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1. Introduction 

This talk summarized some recent work [1] in collaboration with Ted Rogers. 

The overall motivation is to understand the evolution of transverse-momentum-dependent 
(TMD) parton densities (etc) especially at the relatively low values of Q that have considerable 
current interest, as can be seen from many talks in this session on spin physics. We are particularly 
concerned with the non-perturbative part of the evolution, since there appear to be inconsistencies 
between the evolution found in fits [2, 3] to higher Q Drell-Yan data, and the slower evolution pre¬ 
ferred (e.g., [4, 5]) by more recent data at lower Q. When TMD functions are Fourier transformed 
into a space of transverse coordinates Bt;, non-perturbative contributions are at large bj. 

Our aim was to to try to preserve good fits to the Drell-Yan data, while also agreeing with the 
lower energy data and satisfying non-perturbative constraints from the theory side. 


2. Review of TMD factorization and the organization of non-perturbative 
information 


We use the following TMD factorization formula for the Drell-Yan cross section differential 
in the lepton pair momentum and lepton angle: 


da 

d‘^qdD. 



dajj{Q,ix^Q) 

da 


J g“iT bv fj/B{^B,bT-,Q^,Q)d^bj, 


( 2 . 1 ) 


valid when qj Q and polarization effects are ignored. The functions f{x,bT,Q^,Q) are the 
Fourier transformed parton densities (pdfs) to bj, with the CSS and jJ. parameters set to and 
Q. The perturbative har'd scattering factor is dd(2,/i i—)■ Q). The TMD pdfs obey an evolution 
equation of the form 


dln/////(v,z?T;6^;6) 

dlnQ 


fQ du 

r{«s{Q))+K{br,Q) = Y{a,{Q))- —YK{as{l^))+K{bT,Hb)- (2.2) 

JBb F 


The strongly universal function K{bj',lJ.) controls evolution of the shape of the TMD functions, 
and its behavior at large bj is the primary concern of our work. In the right-most part of (2.2), a 
renormalization-group transformation to scale I/b^, was applied to remove large logarithms. 

Non-perturbative information is (a) in the values of the TMD pdfs and of K at large bj, and 
(b) from ordinary pdfs that appear in the OPE that gives the TMD pdfs at small bj. 

To separate non-perturbative contributions to evolution, we use the CSS method to write: 

dln//m(v,(7T;2^;G) fQ du 

- ^ ^ -= 7(«^(2))-/ —(2.3) 

dlnQ ft 


where a smooth cutoff on the perturbative part is provided by Z;* = Z^x/y 1 +^T/^max- Perturba¬ 
tive calculations give the first three terms on the right of (2.3), while fits to data are made for a 
parameterized form for gK{bY,bmay.), which includes the non-perturbative large Z^t contributions. 

The predictive power of this TMD factorization formalism beyond the calculable perturbative 
contributions is from two sources. First is the universality of pdfs between reactions. Second is the 
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strong universality of K, i.e., its lack of dependence on the reaction, on the hadron and quark flavor, 
and on spin and v. The Q dependence of K is governed by the perturbative function Yk, leaving 
non-perturbative dependence as a function of bj only. 

Common choices of are 0.5 GeV^^ and 1.5GeV^^ Observe that if bmax is chosen to be 
too conservatively small, then fitting of gK to data includes reproducing the full in a region 

of bj that is still accessible to perturbative calculations. 


3. Geography of evolution 


E605 Data 




CDF Z Run 1 




Q: 7-18GeV, v^ = 38.8GeV 


e = ISOOGeV 


Figure 1: Left, top; cross-section for Q from 7GeV^* to 18GeV^^ in the E605 experiment compared with 
fits by Landry et al. (BLNY), adapted from [2]. Left, bottom: corresponding bj integrand at g = 11 GeV^* 
for fits by Konychev and Nadolsky (KN) [3]. Right: similar plots for Z production at the Tevatron; note that 
the cross section is for da / dqi instead of da / which gives a kinematic zero at qi = 0. 


As illustrated in Fig. 1 , evolution to higher Q shifts the dominant region of bj to ever lower 
values. Therefore the differential cross section as a function of qi broadens as Q increases. The in¬ 
creasing suppression of the large region implies that the cross section is eventually dominated by 
perturbative effects, even at <77 = 0, provided that the large distance properties are non-pathological. 
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4. Results for K from fits 

The evolution of the integrand W in Eq. (2.1) is given by 
dlnW fQ da 

—^ = G(a,(e)) - / — rKioc,{n))+K{b^-Hb,)-gKibr,brmLx), (4.1) 

ainL; M 

where G{as{Q)) is perturbatively calculable. The right-hand side is the sum of a ^-dependent 
term and a Z^T-dependent term. The ^-dependent term affects only the normalization of the cross 
section. 

The change in shape of the cross section is governed only by the Z^T-dependent part, which can 
be considered as the value of ^(Zjt; Mi) at some fixed reference scale ML 

/■Ml dn 

K{b'Y\ll\) = — j - 7/r(05.s(M))(4.2) 

Jiibt M 

Therefore we can gain an understanding of the change of shape of the ^j-dependent cross section 
with Q from the functional dependence of K{bT', Mi) on bj. 

k(bT\ Q = 2GeV) 



Q: mz lOGeV 2GeV 

Figure 2: K and its components from the KN and BLNY fits. See the text for details. 

In Fig. 2 is shown the Z^t dependence of .^(Z>T;2GeV^^) in the KN and BLNY fits. Overall, 
we see a decreasing function, which corresponds to the shift to smaller dominant values of bj with 
increasing Q that we saw in Fig. 1. 

The red curve gives the purely perturbative, renormalization-group-improved prediction for K. 
This is obtained by setting Zj^ax = °° and g/f = 0 on the right-hand side of (4.2). The calculation 
shown was made with 2-loop approximations for the evolution of the running coupling and for Jk- 
This prediction is valid when bj is small enough, but then diverges to —oo from the Landau pole in 
the approximated coupling, where perturbative-based calculations are completely untrustworthy. 

CSS’s Zj*- prescription gives a smooth cut off of the perturbative term, intended to restrict it 
to a region where perturbatively based calculations are valid. The cut-off part of K, i.e., the first 
two terms on the right of (4.2) is shown as the dotted curves in Fig. 2. The blue curve is for 
^max = 0.5GeV^*, corresponding to the BLNY fit of [2]. The black curve is for Z^max = 1.5GeV^', 
corresponding to the KN fit [3], which gives a better fit to the Drell-Yan data. 

Finally we include the fitted functions gK for the two fits giving the black and blue solid curves. 
In both cases gK is purely quadratic: gx bj. We first notice that each of these two curves matches 
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the purely perturbative red curve beyond where the bmax cut off is important. Indeed the KN curve 
gives a good match to above 2GeV^^ = 0.4 fm, which is, a priori, a region of at best marginal 
applicability of perturbative methods. That the fits match a perturbative calculation suggests that 
the primary result of fitting the one parameter in each fit is to reproduce perturbation theory and 
then to extrapolate the result to large b^: in a non-singular fashion. 

At large there is a dramatic difference between the KN and BLNY curves. This is reflected 
in the factor of two difference between the corresponding integrands at the right-hand edge of the 
lower right plot in Fig. 1. But that large fractional difference is in a place where the integrands are 
small, and so has little effect on the quality of the fits. 

To understand where the differences between the curves matter, we need to know the typical 
values of b'l involved. These are given below the graph in Fig. 2, and we deduce that the fits were 
primarily sensitive to Z^x below about 2GeV^'. The values at larger b^ are only an extrapolation, 
which need not be correct. To probe larger values of Z^x experimentally, we need lower Q. 

That the extrapolation is actually wrong is indicated phenomenologically by the results of Sun 
and Yuan [4]. At large-Z^x, the fitted Gaussian for the TMD pdfs combined with evolution governed 
by a quadratic gK gives a ^-dependent Gaussian behavior for the integrand IT: 


W ^ g-^rkoeffW+constinfeVGo)] = _ _ (at large Z^x). 


(4.3) 


At a value of Q appropriate for data from HERMES and CEAS, the value of a{Q,x) in this equation 
becomes negative when the BENY fit is used; this gives a completely unphysical cross section. The 
KN fit (with Z^max = 1.5GeV^' = 0.3fm) merely gives a value of a{Q,x) much too low to agree 
with the data. 

5. Improved large-Z^x properties 

Phenomenologically, we have good standard fits to Drell-Yan data that determine K for Z^x 
up to around 1.5GeV^^ = 0.3 fm, but not much further. At lower Q, larger Z^x dominates, but 
the extrapolated evolution appears to be too strong to agree with data. Eurthermore, from the 
theoretical side the Gaussian large-Z^x behavior of TMD functions is disfavored [6]. Instead, the 
expected behavior of Euclidean correlation functions is an exponential times a power: 


(5.1) 



That is a non-perturbative statement, with m being the mass of a relevant state. (One could, of 
course, ask whether this expectation is really correct in a confining fheory like QCD) 

Supposing fhaf fhe mass in fhe exponenf is energy-independenf suggesfs fhaf K{bT') goes fo 
a consfanf as Z^x —^ The consfanf gives a Q-dependenf change (presumably a decrease) fo fhe 
normalization of TMD pdfs af large bj, buf does nof affecl fhe exponential ilself. 

We Iherefore proposed [1] a new parameferizafion for gK- 



(5.2) 
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where 


is the value of gx at bj 


2Cr- /'O/t’max dU^ 
go(^max) = g0(Vax,0) + — / 

^ -'Cl/6„ax.O M 

= oo. This is arranged to have the following properties: 


(5.3) 


• At moderate bj it is approximately quadratic, and the sum of gx and the cut-off K{b^',jJLh^) 
approximately agrees with perturbation theory. 

• At large bj, it goes to a constant. 

• The constant is given ^max dependence to compensate (to leading order), the (^max dependence 
of the perturbative term: \mib^^ooK{b^,-,fLb^) = K{bmax',l^b^.^J- 

• Only one free parameter is used. 


As we well see, the result gives much reduced (j^ax dependence of evolution compared with 
standard parameterizations. The exact is a quantity in full QCD and therefore does not itself have 
any b^^x dependence whatsoever. If a simple quadratic function is used and happens to be correct 
for gx at one value of byaax, it cannot be valid for other values of b^ax', a different functional form 
is needed. 

Of course, we do not imagine that our proposed formula (5.2) is exactly correct; we simply 
propose it as a starting guess from which further refinement is possible. 

k(bf. Q = 2 GeV) 



Figure 3: Results of new parameterization for K. 

We have not yet used our new parameterization to perform comparisons or fits with data. To 
show that it is likely to give reasonable results, we show in Fig. 3 what it gives for K for two values 
of bmax, with the choice that go = 0-3 when b^^x = 1.5GeV^^ It gives reasonable agreement with 
the KN results in the region of where the KN fit was determined by Drell-Yan data. But the 
flattening of the curves at larger bj shows that the evolution of the shape of the distribution is 
slower at large Q, as is necessary to be compatible with data at lower energy. 
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From P. Nadolsky, T. Rogers 
and B. Wang 


Figure 4: W with new parameterization, compared with KN version. 

Finally, to measure the compatibility with the Drell-Yan data, Fig. 4 shows some results (from 
Nadolsky, Rogers and Wang) for the Fourier transform in (2.1), atQ = 20GeV, in comparison with 
the KN results for two values of bmax- 

6 . Summary of results 

We argue that K goes to a constant at large bj. This value should be measured, of course. 
We propose a new parameterization to interpolate between the constant at large bj and the known 
behavior at moderate bj- It should give better agreement with data and general principles over a 
wide range of Q. Further refinement is possible, of course. 

To facilitate comparison between different work, fits should be presented in terms of the full 
^{b'v), not just in terms of gK- 

In [1], we argued that the following function can give useful diagnostics: 

( 6 . 1 ) 

It controls evolution of shape of TMD functions, is scheme and scale independent, and is strongly 
universal. 
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